The Cascade of Circulations in Fluid Turbulence 
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Kelvin's Theorem on conservation of circulations is an essential ingredient of G. I. Taylor's the- 
ory of turbulent energy dissipation by the process of vortex-line stretching. In previous work, we 
have proposed a nonlinear mechanism for the breakdown of Kelvin's Theorem in ideal turbulence 
at infinite Reynolds number. We develop here a detailed physical theory of this "cascade of circu- 
lations". Our analysis is based upon an effective equation for large-scale "coarse-grained" velocity, 
which contains a turbulent-induced "vortex-force" that can violate Kelvin's Theorem. We show that 
singularities of sufficient strength, which are observed to exist in turbulent flow, can lead to non- 
vanishing dissipation of circulation for an arbitrarily small filtering length in the effective equations. 
This result is an analogue for circulation of Onsager's theorem on energy dissipation for singular 
Euler solutions. The physical mechanism of the breakdown of Kelvin's Theorem is diffusion of lines 
of large-scale vorticity out of the advected loop. This phenomenon can be viewed as a classical 
analogue of the Josephson- Anderson phase-slip phenomenon in superfluids due to quantized vortex 
lines. We show that the circulation cascade is local in scale and use this locality to develop concrete 
expressions for the turbulent vortex-force by a multi-scale gradient-expansion. We discuss implica- 
tions for Taylor's theory of turbulent dissipation and we point out some related cascade phenomena, 
in particular for magnetic-flux in magnetohydrodynamic (MHD) turbulence. 

PACS numbers: 47.27.Ak,47.27. Jv,47.32.Cc,47.37.-l-q 



I. INTRODUCTION 

The fundamental laws of vortex motion for incompress- 
ible inviscid fluids in three space dimensions were formu- 
lated by Helmholtz |lj . Starting from the incompressible 
Euler equations for an ideal fluid, he showed that vor- 
tex lines are material lines and that the flux within any 
vortex tube is a Lagrangian invariant. Lord Kelvin P] 
gave an elegant alternative formulation of these laws in 
terms of the conservation of circulation, for any closed 
loop advected by an ideal fluid. This theorem is equally 
valid in any space dimension. 

However, all of these results depend upon an implicit 
assumption that the solutions of the fluid equations re- 
main smooth in the inviscid limit. In this limit, as the 
Reynolds number tends to inflnity, all smooth, laminar 
solutions of the Euler equations are unstable and the 
fluid motion becomes turbulent. For inflnite- Reynolds- 
number turbulent solutions, standard conservation laws 
of the ideal Euler equations of motion need not hold. For 
example, both experiments d, ^ IE Hi and simulations 
0,13 show that energy is not conserved in turbulent fluids 
even in the limit as molecular viscosity tends to zero. The 
anomalous rate of energy dissipation in turbulent fluids 
was attributed by Onsager 9] to predicted Holder singu- 
larities in the solutions of the inviscid Euler equations. In 
particular, he showed that a (spatially-minimum) Holder 
exponent /imin < 1/3 is necessary for an Euler solution to 
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dissipate energy. See also 0, 0, 0| . The existence 
of such near-singularities for turbulent velocity fields at 
high Reynolds number has been confirmed by data from 
experiments and simulations 

In a previous work flT^ (hereafter referred to as "I" ) we 
considered similar questions for the conservation of circu- 
lations by turbulent solutions. In that paper we proved 
an analogue of Onsager's theorem, stating necessary con- 
ditions for the anomalous dissipation of circulations by 
inviscid Euler solutions. Furthermore, since these condi- 
tions are expected to be satisfied in turbulent flow, we 
conjectured that Kelvin's Theorem, in its usual form, in- 
deed breaks down for the relevant high-Reynolds num- 
ber solutions. We termed this phenomenon a "cascade 
of circulations." In a following paper we presented 
evidence from direct numerical simulations for the ex- 
istence of such a cascade. The purpose of the present 
paper is to elaborate further the physical theory of this 
phenomenon. In particular, our aims are as follows: 

In the remainder of this section of the paper, we shall 
discuss some important background information. We 
first remind the reader of the classical Kelvin Theo- 
rem. Next we briefly review some ideas of G. I. Taylor 
[l^ I20I [2H about the role of circulation-conservation in 
the production of energy dissipation in three-dimensional 
turbulence. In the second section of the paper we present 
our new results. First, we discuss the flltering approach 
which is the basis of our theory, and explain its relation to 
renormalization-group (RG) ideas and to large-eddy sim- 
ulation (LES) modelling of turbulent flows. Second, we 
establish exact results for large-scale circulation balance 
of low-pass filtered velocity fields. Third, we explain how 
Taylor's argument can be extended to stretching of fil- 
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tered vorticity and how this is related to forward cascade 
of energy through the inertial range. Fourth, we review 
the resuhs from I on the possibihty of anomalous dissi- 
pation of circulations in the limit of zero filtering length. 
Fifth, we point out an interesting analogy between this 
cascade of circulations and the phenomenon of phase-slip 
in superfluids, noting similarities with previous ideas of P. 
W. Anderson ^22] . Sixth, we discuss the scale-locality of 
the circulation-cascade and elaborate a multi-scale gradi- 
ent (MSG) expansion for circulation-flux, along the lines 
laid out earlier for turbulent stress [23L 12^ . Finally, in 
the conclusion section we discuss some implications of our 
results and various extensions to magnetohydrodynamic 
(MHD) and geophysical fluid turbulence. 

A. Classical Kelvin Theorem 

We here briefly review some standard facts about the 
conservation of circulations. Let u(x,t) be a smooth 
velocity field solving the incompressible Navier-Stokes 
equation with viscosity v 

^fU -I- (u-V)u = - Vp + A u, V-u = (1) 

where x G A C R'^, for any integer d > 2. Here p(x, t) 
is the so-called "kinematic pressure" (or, thermodynam- 
ically, the enthalpy per unit mass). For any closed, ori- 
ented, rectifiable loop C C A at an initial time to, one 
defines the circulation 

r{C,t)^(f) u{t)-dx^ u:{t)-dA (2) 

Jcit) Js{t) 

where C{t) is the loop at time t advected by the fluid 
velocity, S{t) is any surface spanning that loop, and 
uj{t) — Vxu(t) is the fluid vorticity. These circula- 
tions satisfy the Kelvin-Hclmholtz theorem in the fol- 
lowing sense: 

^T{C,t)^vi Au(t)-dx. (3) 

dt Jc(t) 

E.g., see [2^, section §1.6, for the standard derivation. It 
is worth observing that the Kelvin theorem for all loo ps C 
is formally equivalent to the Navier-Stokes equation |26l |. 
Indeed, if u(x, t) is a smooth spacetime velocity field, 
divergence-free at all times i, then equation (jJl implies 
that 

<j> [Dtn{t) -i^A u{t)] -dx = (4) 

for all loops C at every time t. Here DtU — 9iU+ (u-V)u 
is the Lagrangian time-derivative and the equation Q is 
derived by applying © to the pre-image of the loop C 
at initial time to. By Stokes theorem, equation can 
hold for all loops C C A if and only if there exists a 
pressure-field p(x, t) such that the Navier-Stokes equa- 
tion holds locally and also globally, if the domain A 
is simply connected. 



In the inviscid limit v 0, the circulation is formally 
conserved for any initial loop C. The fluid equations in 
this limit, the incompressible Euler equations, are the 
equations of motion of a classical Hamiltonian system. 
They can be derived by the Hamilton-Maupertuis prin- 
ciple from the action functional 

S'[x] = i [ ' dt [ dsL \±{a,t)\^ (5) 

2 J to J A 

with the pressure fleld p(x, t) a Lagrange multiplier to 
enforce the incompressibility constraint. Here x(a, t) 
is the Lagrangian flow map which satisfies x(a, t) = 
u(x(a, t\ t) with initial condition x(a, io) = a- See 
|27l |28| for reviews. This variational principle yields 
the fiuid equations in a Lagrangian formulation, as 
x(a, i) = — Vp(x(a, t), t). The Eulerian formulation 
(with 1/ = 0) is obtained by performing variations in 
the inverse map a(x, t), or "back-to-labels map", with 
fixed particle positions x. This Hamiltonian system has 
an infinite-dimensional gauge symmetry group consist- 
ing of all volume-preserving diffeomorphisms of A, which 
corresponds to all smooth choices of initial fiuid particle 
labels. In this framework, the conservation of the circu- 
lations for all closed loops C emerges as a consequence of 
Noether's theorem for the particle-relabelling symmetry 
|29|. For reviews, see Section 4 or [2i|, Section 2.2. 



B. Circulation and Turbulent Energy Dissipation 

In several papers 0, 0, |^, G. I. Taylor has ar- 
gued for the importance of conservation of circulations 
in the turbulent generation of energy dissipation at high 
Reynolds numbers in space dimension d — 3. We briefiy 
review his ideas. The simplest version of Taylor's argu- 
ment is based upon the concept of vortex line-stretching. 
Consider a vortex tube initially with length Lq, cross- 
sectional area Ao, and vortex strength wq. Taylor as- 
sumed that such a vortex tube at high Reynolds num- 
ber will evolve as a material line. Taylor also reasoned 
that vortex lines (or any material lines) should tend to 
lengthen, on average, under random advection by a tur- 
bulent velocity field. Thus, at a later time t > to, the 
tube length is typically L{t) > Lq- By incompressibil- 
ity, the volume V{t) = L{t)A{t) does not change in 
time, so that A{t) < Aq. Furthermore, Taylor reasoned 
by the Helmholtz theorem that the vortex-fiux through 
the tube, T{t) = uj{t)A{t), would not change, so that 
Lu{t) > uJo- In fact, by this chain of reasoning, 

Lo{t)/ujo = L{t)/Lo (6) 

and vortex strength increases in direct proportion to line- 
length. Because the viscous energy dissipation in the 
vortex-tube is given by / w'^{t)dV — vuj'^{t)Vo, this 
process should lead to a dramatic enhancement of dis- 
sipation. 
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However, this argument contains an apparent incon- 
sistency. On the one hand, Taylor's assumptions that 
vortex-hnes are material lines and that the Kelvin Theo- 
rem applies require that the viscosity term in the circu- 
lation balance ^ can be neglected. On the other hand, 
Taylor retains the viscous dissipation in the energy bal- 
ance, arguing, in fact, that it is sizable. It is not at all 
clear that it is valid to ignore the viscosity effects in one 
place and to keep them in another. Taylor himself recog- 
nized the delicacy of his argument. In 20] he presented 
this line of reasoning, and then wrote: "When oj^ has 
increased to some value which depends on the viscosity, 
it is no longer possible to neglect the effect of viscos- 
ity in the equation for the conservation of circulation, so 
that (10) [our ijHJl] ceases to be true." Thus, Taylor as- 
sumed that there is some interval of time or some range 
of length-scales for which viscous effects can be neglected 
in the circulation balance Q. We shall critically review 
this assumption below. 

In a following paper l21j, Taylor tested some predic- 
tions of his argument using experimental data for decay- 
ing turbulence generated from a wind-tunnel. His analy- 
sis was based upon the following equation for production 
of enstrophy. 



v\Vu}\'^ 



(7) 



which is an exact consequence of the incompressible 
Navier-Stokes dynamics ^ for space-dimension d = 3. 
Here Sij = {l/2){dui/dxj + duj/dxi) is the strain ma- 
trix. Under conditions of space-homogeneity, the average 
of the transport term vanishes, so that 



(dM)(i|a;|2) 



(8) 



Taylor's argument on vortex-line stretching suggests that 
{d/dt){\uj\^) > 0, which can hold if and only if {iv^Siv) > 
i^(|Va;|^) > Thus, enstrophy will be created when the 
mean rate of vortex-stretching by the strain is positive 
and exceeds the mean destruction of enstrophy by viscos- 
ity. In [23| , Taylor found from an analysis of wind-tunnel 
data that the latter condition holds for an initial range 
of time in decaying turbulence. 



II. CIRCULATION CASCADE 

We now turn to an analysis of circulation conserva- 
tion in high-Reynolds-number turbulent flow. One ap- 
proach would be to directly analyze the z/ — > limit of 
eq. (|2Jl. However, we shall pursue a complementary ap- 
proach based upon a study of nonlinear transfer in the 
inertial range. 



A. Filtering Approach 

To analyze the dynamics in the inertial range, we in- 
troduce effective equations that govern the evolution of 
the velocity field at large length-scales. For any chosen 
length £, let 



u^(x) = J drGe{r)u{x + r) 



(9) 



denote the low-pass filtered velocity at scale £, where 
Ge{r) = i~'^G{r/£) is a filter kernel. We shall assume 
that G is positive, smooth, rapidly decaying in space and 
with unit integral. Then ug satisfies an effective equation: 



dtue + (urV)u£ + V-r^ = - Vp^ + iy Aue, (10) 

where is the filtered pressure and is the turbulent 
stress-tensor 



Te = (uu)^ - UiUi. 



(11) 



The filtering operation that we have employed can 
be regarded as a "coarse-graining" that eliminates 
high-wavcnumber modes, as in renormalization-group 
methodology |^ . Because of momentum conserva- 
tion, the effective renormalized equation can change only 
by additional contributions to the stress tensor. This 
filtering approach is also the mathematical basis of the 
large-eddy-simulation (LES) modeling scheme [13, Issf . 
In this scheme, the stress tensor is the main unknown 
which must be modelled, in order to obtain a closed equa- 
tion for computation of the large-scale velocity field. 

In the inertial-range of turbulent flow the final viscosity 
term in ea. ljlUII can be neglected. For example, a fairly 
crude estimate based upon the identity 



Au£(x) 



dr(AG)f(r)u(x + r) (12) 



is A ue\\2 < (j^/£^)(const.)||u||2, where ||u||2 = 
[/ dx|u(x)p]^/^ is the L^-norm. If the total kinetic en- 
ergy per mass E = (l/2)||u||| remains finite in the limit 
as — > 0, then the viscosity term in ea. pO|l tends to zero 
in L^-norm for any fixed filter-length £. 

There is another form of the effective equation l|10(l 
which is useful. Note that the stress appears only via the 
turbulent (suhgrid) force ii — —'V-ti [33 ]. This can be 
replaced in (|10|) using the following elementary identity 



ff. = -Vkf 



(13) 



where kg = (l/2)trTf is the turbulent kinetic energy |3S 
and 



(14) 



is the turbulent vortex force. With this replacement, H10() 
becomes 



dtu-i + (u£-V)u£ = -Vpl + fg* + A U£, 



(15) 
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where =Pi + ki is a modified pressure. Although this 
form of the large-scale effective equation leads to more 
intuitive results, it is less easy to make sense of math- 
ematically. In fact, the vortex force f^* could be badly 
ultraviolet divergent in the limit as ^ 0. Notice that 
for infinite-Reynolds-number turbulence the velocity u is 
believed to be a continuous but non-differentiable func- 
tion, so that the vorticity u: exists only as a distribution. 
Therefore, the product uXu; is a priori ill-defined. How- 
ever, the vortex force remains well-defined due to the 
identity H13|) . since both and Vfc^ make sense as long 
as ||u||2 < oo. 

B. Circulation-Balance in the Large-Scales 

It is natural to inquire about the circulation-balance 
for the large-scale effective equation. Let us choose an 
oriented, rectifiable, closed loop C in space. We define 
Ci{t) as the loop C advected by the filtered velocity u^. 
This definition makes sense, since the filtered velocity 
ue is Lipschitz in space, and corresponding flow maps 
Xf(a, i) defined by 

{d/dt)xi{a,t) ^Ui{xi{a.,t),t), x^(a,io)=a, (16) 

both exist and are unique (see I) . We define a "large-scale 
circulation" with initial loop C as the line-integral 

TeiC,t)^<f_ Mt)-dx. (17) 

Jctit) 

for i < R = the radius of gyration of the loop C [tJI . The 
same calculation that establishes the Kelvin theorem, but 
using the effective eq. ()inil rather than Navier-Stokes 
eq-O, gives 

{d/dt)Te{C, t)^ i [h{t) + A uf] • dx (18) 
Jcdt) 

If the Navier-Stokes eq.(^ were driven by an external 
body- force P'^', then there would be an additional term 
f inside the square bracket in ea. H18|l . If this external 
force is spectrally supported at wavenumbers of order 
1/L, then its contribution to the circulation balance is 
0{R/L). Thus, the forcing term is neghgible for R <^ L. 
Likewise, the viscous term in ea. (|18|l is negligible for 
small viscosity ly and fixed filter-length £, by an elabora- 
tion of the argument given around ea. lfT^ . (A so-called 
"trace theorem" can be used to estimate the restriction 
of Aue to the loop C; see I and H^l). 

These remarks show that the nonlinear term from the 
subgrid force is the dominant term in the circulation bal- 
ance (I18|l for incrtial-range values L R > £ ^ rjd 
(where 77^ is a dissipation length-scale determined by the 
viscosity . If we imagine that the total circulation at all 
scales on the loop is conserved, then the line-integral of 
on the RHS of H18|) represents a "transfer" of circulation 
to subgrid modes at length-scales < £. This motivates 



the definition, for any loop C and filter length £, of a flux 
of circulation 

Ke{C, t) = - (f_ f^(t) ■dx = - (f_ f;{t) ■ dx (19) 

JCi(t) Jci{t) 

so that {d/dt)Ti{C,t) = —Kii{C,t) (up to small correc- 
tions from external forcing and viscosity) . We have used 
identity H13I) to justify the equality of the two expres- 
sions m the definition UH). The minus sign has been 
introduced so that the signs of the circulation H17(l and 
the circulation-flux H19() should be positively correlated. 
This expectation will be discussed more below. 

The "circulation-flux" defined in H19() has the physical 
dimensions of work or of torque (per unit mass) . Addi- 
tional insight into its meaning can be obtained by decom- 
posing the turbulent vortex-force l|14|l into components 
perpendicular and parallel to large-scale vortex-lines: 

= cTtXCji, f||^ = (f/-d>£)a)£ (20) 
where Cj^ ~ u3i/\'Ujt\ and 

= c2;,xf;. (21) 

If ti is the unit tangent vector to the curve Ci{t) and s 
is the arc- length parameter, then 

Ki{C,t) ^ (f)_ (Ti{t)-nids - (f)_ f|^(t)-dx (22) 

where = t^Xw^. Note that the latter vector is normal 
both to lines of large-scale vorticity ZU^ and to the loop 
Ci{t), but it is not generally a unit vector. The first term 
in H22|) can be interpreted as a lateral diffusion of vortex- 
lines out of the advected loop, where erg plays the role 
of a transport vector of vortex-lines. The second term in 
(I22|l represents an additional work (or torque) due to the 
parallel component of the turbulent vortex-force. 

Some particular cases of (|^ are of special interest. 
For example, consider the case that Ce{t) is instanta- 
neously a closed vortex line. (This property will not gen- 
erally be preserved in time). Then the first term in H22() 
vanishes and Ki{C,t) = — ^^^^^ f|*^(t) • dx. Such inte- 
grals play an important role in vortex-reconnection the- 
ory l35j. The distinguished vortex lines for which this 
integral is extremal drive the reconnection process and 
the value of the integral for such lines gives the rate of 
reconnection of vortex-flux. This integral is therefore the 
proper point of departure for a theory of turbulent re- 
connection of large-scale vortex-lines. Another special 
case of (|22|) of interest is when the loop Ci{t) lies in a 
transversal surface normal to the lines of large-scale vor- 
ticity. In that case, the second term in (|22|) vanishes and 
Ke{C,t) — f^^^^-^(Te{t)-hids, where hi — igXuJi is now 
a unit vector. This condition is always satisfied for space 
dimension d = 2. The flux of circulation is then entirely 
due to the diffusion of vortex-lines out of the loop. 

These remarks on physical interpretation of Ki{C,t) 
lead to some natural guesses on the correlation of its 
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sign with that of the circulation Ti{C, t). The latter can 
be written as 



uJe{t)-dA, 



(23) 



Se{t) 



where Si{t) is any smooth surface spanning the loop 
Ce{t) and with orientation consistent to that of Ci{t) 
(by the righthand rule). If the circulation is positive, 
then there is a net contribution from vortex lines thread- 
ing the loop in the direction of the surface unit normal. 
If the effect of the turbulence is "diffusive" on average, 
then one would expect that the vortex-force will tend to 
smooth out the excess of positive-sign vorticity thread- 
ing the loop. Thus, according to the sign convention of 
the definition H19|l . we can expect that Ki{C, t) will also 
tend to be positive and to reduce the overall magnitude 
of the large-scale circulation. Of course, this argument 
works equally well when Vi{C,t) has negative sign. We 
may therefore expect that there is in general a "forward 
cascade" of circulations, and that the magnitude of the 
large-scale circulation, of whatever sign, will tend to be 
decreased by the small-scale turbulence. This reasonable 
result has been confirmed by numerical results in [T^ . 

An interesting exception is the inverse-energy cascade 
for d = 2 turbulence. For space dimension d ~ 2, the 
enstrophy il(t) = {1/2)\lu\'^ is an inviscid invariant and 
its flux to unresolved scales < is measured by 



(24) 



where oj^ is the filtered vorticity (perpendicular to the 
plane) and er e is the vorticity transport vector defined in 
(EH- See |3ll37ll3^. From ^ one can see that enstro- 
phy will cascade forward to small scales when vorticity 
transport tends to be "down-gradient" and VcU^ -erf < 0. 
On the other hand, enstrophy flux will be inverse to large- 
scales when the vorticity transport is "up-gradient." In 
d = 2 there are expected to be two inertial cascade 
ranges, the direct enstrophy cascade where the mean 
enstrophy flux is positive and the inverse energy cas- 
cade where the mean energy flux is negative |39l l4Clj| . 
However, there is also some "leakage" of energy flux 
and enstrophy flux into the opposite ranges (e.g. see 
HE 01)- In particular, the mean enstrophy flux in the 
inverse energy cascade range is negative^ or toward larger 
scales. This means, according to (|24|l . that the vorticity 
transport in that range is, on average, "up-gradient" or 
"anti-diffusive" . Therefore, our argument for the sign 
of circulation-flux is reversed. In the inverse cascade 
range, a loop containing an excess of one sign of vortic- 
ity should tend to accumulate more vorticity of the same 
sign. Thus, in the d = 2 inverse energy cascade range 
there should be also an "inverse cascade of circulations" 
III. 



C. Stretching of Large-Scale Vorticity 

We have seen that the "large-scale circulations" , in the 
inertial range, evolve according to the equation 



{d/dt)Ti{C,t) = i r^{t)-dyi. 

JCiit) 



(25) 



The term on the righthand side due to the vortex-force 
need not be negligible. Thus, Taylor's conjecture that 
Kelvin's theorem should hold in the inertial range, even 
approximately, is far from obviously true. In the next 
section we shall explore this question mathematically, to 
the extent possible. Here we discuss some physical im- 
plications of Taylor's conjecture, if true. 

If we suppose that the inertial-range circulations are 
conserved, then Taylor's argument about vortex line- 
stretching can be repeated for filtered vorticity, implying 



id/dt){\LJe\'') >0. 



(26) 



This result can also be understood from the equation 
for the filtered vorticity, obtained by taking the curl of 
equation H15|) (with ly = 0): 



dtujt, + {ue-V)uJe = {uJe-V)ue + Vxf/ 



(27) 



From this an equation for inertial-range enstrophy easily 
follows: 



1 



^ 2 



1 



a;; SeLJi + • ( V X tJ^) . (28) 



[Compare with |4j|, eq.(51), for ly — > 0.] Notice that the 
vorticity transport vector ere defined in contributes 
to the space transport of enstrophy. However, assuming 
space-homogeneity, all of the space-flux terms average to 
zero and 



(d/dt) {^\ZJef) = {ZjJSeZJe) 



(f;-(VxaJ,)). (29) 



This equation is an exact inertial-range analogue of equa- 
tion (jHl for total enstrophy. The first term on the 
righthand side of (|29() represents inertial-range vortex- 
stretching and the second term represents enstrophy 
flux to length-scales < £. For freely decaying turbu- 
lence at early times, Taylor's argument predicts that 
(ZU|S,ZUf) + (f;.(VxZU,)') = id/dt){^pe\^) > 0. On 
physical grounds, one expects that the vortex-stretching 
is positive and the enstrophy transfer term negative, 
with the net enstrophy production positive. At later 
times a quasi-equilibrium should be established so that 
{d/dt){^\ZJe\'^) « and the dominant balance becomes 



< {ujJSeZJi) « -(f;-(VxaJ£)) 



(30) 



For some experimental results on these questions, see 

m. 
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It was observed in |43 that the energy flux He to un- 
resolved scales < i can be expressed approximately in 
terms of the negative skewness of filtered strain and the 
stretching rate of filtered vorticity: 



-tr 



(S,') + 



(31) 



This expression is the first term in a systematic "multi- 
scale gradient expansion" . It follows from an identity 
of Betchov that for any homogeneous turbulence 



(32) 



Thus, the energy cascade will be forward to small scales 
when the mean rate of vortex-stretching is positive. This 
is an inertial-range version of Taylor's mechanism [20l 



D. Anomalous Conservation of Circulation 

We now consider the question whether Kelvin's The- 
orem can hold, in any sense, in turbulent flow at high 
Reynolds number. In view of equation (|18|l or (|25|l . we 
must estimate the magnitude of the circulation-flux de- 
fined in H19|l. The following simple identity, observed 
in ^3 1 is useful to provide an estimate of the turbulent 
subgrid force: 



fiti^) ^ J I dr {djG)fXr) Su,{r;x)6uj{r;x) (33) 



dr{djG)iir)6u,{r;x) / dr' Geir') Suj{r';x) 



Here 6u{r; x) u(x + r) — u(x) is the velocity-increment 
with separation vector r at location x. An upper bound 
easily follows that |ff| = 0{\Su{£)\'^ /£), where Su{i) is 
the maximum magnitude of the velocity-increment for 
separation vectors with |r| < ^ |T7l |. 

If the velocity field were smooth, then |<5m(^)| ^ 
(const. )£ for small £ and the subscale force would van- 
ish as £ ^ 0. However, a turbulent velocity field does 
not remain smooth in the limit as the Reynolds number 
tends to infinity. Instead, theory, simulations, and ex- 
periment indicate that the velocity field is only Holder 
continuous with exponent < h < 1 : 



|(5u(r;x)| = 0{r^'). 



(34) 



At each point x one refers to the maximal value h for 
which (|34|l holds as the Holder exponent at that point. 
There is a spectrum of such singularities in the flow, 
with exponent h occurring on a set S{h) with fractal 
dimension D{h). It was pointed out by Onsager that 
the smallest exponent /imin must be < 1/3 to explain 
non- vanishing energy dissipation in the inviscid limit. 
Parisi and Frisch |45l | invoked a multifractal spectrum 
D{h) of singularities to explain the anomalous scaling of 
pth moments of velocity-increments (so-called pth-order 



structure- functions). Such multifractal spectra of Holder 
exponents have been confirmed by analysis of data from 
experiments and simulations [lj,[il,[i3- Of course, at fi- 
nite Reynolds numbers there are only "near-singularities" 
in the inertial-range of scales and the velocity is smooth 
in the dissipation range, where effects of viscosity are 
important. 

From our estimate below ea. ()33|l . we see that jf^j — 
Qf^£2h-i-^ at any point with local Holder exponent h. 
Thus, the circulation fiux Ke{C, t) will go to zero as£ 
if the smallest velocity Holder exponent /imm is > 1/2 and 
if also the curve C{t) has finite length [131 ■ This is an 
exact analogue for circulation flux of Onsager's result 
for vanishing of energy fiux when h^in > 1 /3. Only a suf- 
ficiently rough velocity field can provide a transport of 
vortex lines which is non- vanishing in the limit as ^ 0. 
However, high Reynolds turbulence in space dimension 
d = 3 has a plethora of singularities with exponents 
h < 1/2. For example, the most probable exponent 
with D{h^) = 3 has a value /i* = 1/3, very close to the 
mean- field Kolmogorov value 0, 0, 0| . Furthermore, 
the curves Ci{t) advected by the large-scale velocity Uf 
are expected to approach a fractal curve C {t) in the limit 
as € ^ Thus, circulation-flux is not likely 

to vanish as the filtering length decreases through the 
inertial-range. Numerical simulations of high-Reynolds- 
number turbulence for d — 3 confirm this prediction 1^. 

There is an important subtlety in the formulation 
of Kelvin's theorem for infinite- Reynolds-number turbu- 
lence that must be mentioned at this point. Recent 
work on an idealized turbulence problem — the Kraich- 
nan model of random advection [i^ — has shown that La- 
grangian particle trajectories x(t), x'(t) can explosively 
separate even when xq = Xg initially, if the advecting ve- 
locit y fi eld is only Holder continuous and not Lipschitz. 
See [43 • Mathematically, this is a consequence of the 
non-uniqueness of solutions to the initial-value problem, 
while, physically, it corresponds to the two-particle tur- 
bulent diffusion of Richardson [s^ • It has been rigorously 
proved in [sil Is^ that there is a random process of La- 
grangian particle paths x(t) in the Kraichnan model for a 
fixed realization of the advecting velocity and a fixed ini- 
tial particle position. This phenomenon has been termed 
spontaneous stochasticity 53] and it is likely that it holds, 
not only in the Kraichnan model, but also for singular 
solutions of the inviscid Euler equations. If so, then the 
advected curves C (t) that appear in the definition of cir- 
culation ^ are likely to be random fractal curves! 

If these speculations are correct, then the time-series 
of circulations r(C, t) are also a stochastic process, for 
a fixed turbulent velocity field. In |17| we have pre- 
sented some plausibility arguments in favor of the fol- 
lowing "martingale property" for this random process of 
circulations: 

{r{G,t)\r{G,T),T <t') = r{c,t'), iort>t'. (35) 

Here (•) denotes the expectation over the ensemble of ran- 
dom Lagrangian paths and we have conditioned on the 
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past circulation history {T{C,t),t <t'}. Heuristically, 

{d/dt){r{c,t)\r{c,T),T <t')^ 

-\ira{Ki,{C,t)\T{C,T),T < t') =0. (36) 

The circulation-flux in is conjectured to average to 
zero, due to increasingly rapid oscillations of the vortex- 
force around the loop Gi{t), as ^ ^ 0. See 0|. The 
result in (|36(l has been partially confirmed by the results 
of a numerical simulation in JJ] , providing some support 
to the conjecture H35|l . This "martingale property" is a 
statement of conservation of circulations, in a conditional 
mean sense. It is not clear yet whether this weakened 
version of the Kelvin theorem is valid and, if so, whether 
it suffices for Taylor's vortex-stretching mechanism. 

E. Analogy with Phase-Slip in Superfluids 

It is worth pointing out an analogy of the "circula- 
tion cascade" discussed above with another physical phe- 
nomenon, the "phase-slip" due to quantized vortex lines 
in superfluids [l^, 0| • Anderson had already discussed 
classical analogues of quantum phase-slip in |22|, Ap- 
pendix B. His starting point was the classical Euler equa- 
tions for an incompressible fluid, written as 

atu = -V/i + uXw, (37) 

where ft, = p+(l/2)|upis the enthalpy. Anderson con- 
sidered the line-integral of the fluid velocity u along a 
stationary curve C connecting two points Pi and P2, 
showing that 

{d/dt) / u(<)-dx = -Ac/i+ / (rfxxu)-u;. (38) 
Jc Jc 

Here A^/i = h{P2) — h{Pi) is the difference of h along 
the curve C. Denoting time-average by (•), this relation 
yields 

A^T^^ / ((ixxu)-a;. (39) 
Jc 

Since vortex lines for smooth solutions of the classical Eu- 
ler equations move with the particle velocity u = d'x./dt, 
the righthand side of (|39|l can be interpreted as an av- 
erage rate of flow of vorticity across the curve C. This 
flow rate is thus equal to the average enthalpy difference 
along the curve. After deriving H39|l . Anderson wrote 
[25|: "We see immediately that this equation is far more 
important in a superfluid, where vorticity is conserved 
and quantized, than it is in ordinary fluids, where in a 
laminar flow, for instance, the right-hand side has little or 
no special significance." One critical difference between 
classical fluids and superfluids is that, in the former, the 
vortex-lines for laminar solutions move with the fluid. 
Thus, if one instead considers a material curve C(t), ad- 
vected by the fluid velocity u, then one obtains 

{d/dt) [ u(i).rfx- A^^^^A (40) 

JC{t) 



with A = (l/2)|up -p, rather than The nontrivial 

term associated to flow of vorticity across the curve is now 
absent and ea. H4U|) for a closed loop yields the classical 
Kelvin Theorem. 

Nevertheless, we have found that it is possible for tur- 
bulent flow to yield a nontrivial result. In fact, by filter- 
ing the Euler equation H37|) one obtains 

dtVLi = -Vhi, + u,,xu}i + f; , (41) 

with the additional vortex-force term. This equation is 
equivalent to 

Dm = dm + (urV)u<> = -Vpl + f;, (42) 

which is our old cg. HlSI) for = 0. As we have seen in 
our earlier discussion of the large-scale circulation bal- 
ance, eq. ()18|l or (|25|l , the turbulent vortex- force provides 
a nontrivial transport of vorticity across material curves. 
Here it is crucial that the velocity field be sufficiently sin- 
gular, to permit a transport which is non-vanishing for 

— > 0. If instead the flow were smooth and laminar, then 
fg* in that limit and filtering the equation would 
lead to no new result. For singular solutions the Euler 
equation H37|) must be filtered to make sense, as a matter 
of principle. In the presence of singularities the equation 
is interpreted in the sense of distributions, which means 
that it must be smeared with smooth test functions. 

Nontrivial results are also possible in superfluids, for 
similar reasons. The superfluid phase order parameter 
^(p^beys the Josephson- Anderson frequency equation |2^ 

hd(p/dt = -(/i + ^mul), (43) 

where is the chemical potential and = (?i/m)V(/5 
is the superfluid velocity. It is straightforward to derive 
from (|43|l the superfluid equation of motion 

DtUs = [dt + Us- V)us = - V(/i/m) - Xu;^. (44) 

Here the final term contains the superfluid vorticity u)s = 
VxUs which is, formally, a delta-function supported on 
singular vortex lines (zeroes of the superfluid density). 
Equation l|44|) is the basis of derivations of Kelvin's the- 
orem for superfluids, e.g. see [s^ in the context of the 
zero-temperature Gross-Pitaevskii equation. Note, how- 
ever, that such derivations require that the advected loop 
not pass through singular points where the superfluid ve- 
locity is ill-defined. Since the quantized vortex lines are 
not material lines in general (e.g. see js^ [s^, El]), it 
is possible for them to mig rate out of an advected loop. 
Examples are given in of the failure of Kelvin's the- 
orem due to the intersection of loops with singularities 
that are, formally, represented by the rightmost term in 
eg . (|44|l . That equation is thus analogous to ea. (|42ll for 
classical turbulence. 

One of the concrete manifestations of quantum phase 
slip is the decay of "persistent" superfluid flow in a thin 
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toroidal ring. E.g. see [59j and references therein. This 
process has a number of similarities to the "cascade of 
circulations" in turbulent flow. The decay of the super- 
flow is mediated by the (thermal or quantum) nucleation 
of quantized vortices which migrate out of the ring. The 
passage of a vortex across the toroidal cross-section in- 
duces by phase-slip a pulse of torque which decreases the 
circulation around the ring. The reduction in the angular 
momentum of the superfluid condensate is balanced by 
a gain in the normal fluid excitations, acting as an an- 
gular momentum reservoir. In the turbulent circulation- 
cascade, the large-scale vortex lines are also not mate- 
rial, because singularities in the velocity fleld allow them 
to diffuse relative to the fluid. The subscale modes at 
length-scales < I act as a reservoir, whose feedback on the 
resolves scales > t provides the vortex-force that drives 
the diffusion. Unlike in superfluids, this is a continuous 
process, since classical vortices are not quantized. There 
is also no need for the singularities to be nucleated as 
fluctuations, since they are everywhere present in the tur- 
bulent flow. Finally, if the "martingale" conjecture H35|l 
is correct, then the turbulent diffusion of vortex-lines is 
not persistent in scale, on average, and does not lead to 
irreversible mean decay of circulations. 



F. Scale-Locality and MSG Expansion 

We have referred to this turbulent diffusion of vorticity 
as a "cascade" of circulations, but we have not shown 
that the process is a local-in-scale cascade. Here we shall 
examine this issue, following the general approach in |23j |. 

We note first that the turbulent vortex-force de- 
fined in (|14() is a priori not ultraviolet (UV)-local, under 
conditions realistic for turbulence in d = 3. In fact, the 
vorticity is a dissipation-range variable and its largest 
contributions come from the viscous scale. The argu- 
ments in [2^ for UV-locality would apply to if the 
Holder exponents hu of velocity and of vorticity both 
were positive. However, /i^^ = /i„ — 1, so that vortic- 
ity is expected to have negative Holder exponents in the 
infinite- Reynolds-number limit (and thus to exist only as 
a distribution) 0. It is possible that there could be 
cancellations in the average (|14|l over displacement vec- 
tors that defines the vortex-force. E.g. this was found to 
be true in the d = 2 enstrophy cascade, by an analysis of 
the results of a numerical simulation [Sg- However, the 
UV-divergence is more severe for d = 3, so that sufficient 
cancellation is less likely there. 

On the other hand, because of the identity (|13|l . we 
may use the turbulent subscale force = — V-r^ rather 
than the vortex-force f^* to study the circulation-flux. 
The force ii has much greater chance to be scale-local, be- 
cause it is defined only in terms of velocity. Indeed, some 
locality properties follow directly from the representation 
(|33|l in terms of velocity increments. As in [23, let us de- 
fine = Ga * u to be the low-pass filtered velocity 
at length-scale A > ^ and define u^"* = u — u^"^ to be 



the high-pass filtered velocity at length-scale 5 < We 
can then define a very large-scale contribution f^-*'^ to the 
turbulent force by replacing both u in the formula H33|) 
with u>^. Likewise, we define a very small-scale contri- 
bution by replacing both u with u'-'^. Now suppose 
that the velocity field has Holder exponent ft, at a consid- 
ered point X. Then, the following estimates can be easily 
derived, by the same methods as in 

|f>^| =0(M2''-2) (45) 

and 

\ff\=0{5'"'/t). (46) 

The estimate (|45ll expresses infrared (IR)-locality. In 
fact, when h < \, this estimate shows that f^'*^ decreases 
for increasing A and fixed I. Relative to the estimate 
|f^| = 0(^2h-i^^ ^Yie estimate ^ for |f>^| is smaller 
by a factor 0((£/A)2(i-'')). Likewise, the estimate (O 
expresses UV-locality. When /i > 0, this estimate shows 
that if-^ decreases for decreasing 5 and fixed I. The es- 
timate (|46|) for If^"^^! is smaller than that for |ff| by a 
factor of 0{{6 / tY^). These results show that most of the 
turbulent subgrid force ig comes, pointwise, from pairs of 
velocity modes at length-scales ~ 

The above arguments do not quite settle the issue of 
locality of the circulation-flux Kg{C,t), however. The 
delicate point here is that large cancellations are expected 
in the line-integral of Eg that deflnes that flux. In order 
to infer scale-locality of Ke{C,t), one must assume that 
similar cancellations occur in the line integrals of f^'**^ 
and f^^. This issue is hard to address mathematically 
but may be investigated using data from simulation or 
experiment. 

The UV-locality properties of the subgrid force ff may 
be used to develop an analytical expression for it, by 
means of a multi-scale gradient expansion jj| . We con- 
sider only the lowest-order term in that expansion, which 
corresponds to the so-called "nonlinear model" for the 
stress | ,33j : 

Ty Cfui^iUj^i. (47) 

Here C = J dr|G(r)prf and a spherically-symmetric fil- 
ter function is assumed, so that ri could be replaced with 
any other single component . (In terms of the constant 
C2 employed in [2^, C = C2/d where d is the space di- 
mension.) We use the convention of subscript ",j" to 
denote dj, so that, for example, u^.j = dui/dxj. We also 
employ the Einstein summation convention for repeated 
indices. To avoid an excess of subscripts, we drop above 
and hereafter the subscript I, since a fixed filter length 
will be always understood. The physical assumption be- 
hind the formula l|T7|) is strong UV-locality, so that only 
adjacent subscale modes contribute to the stress. We ex- 
pect that this extreme assumption is fairly good in the 
d = 3 energy cascade and the d = 2 direct enstrophy cas- 
cade. However, we present arguments below that it fails 
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badly for the d = 2 inverse energy cascade. Note that it 
is already known that the energy transfer is only weakly 
scale-local in d = 2 60iili3. 

From the formula H47|l for the stress, one obtains the 
corresponding formula for the subscale turbulent force: 



(48) 



By means of a standard vector calculus identity, this can 
be written for d < 3 as: 



= C£ tijkUj^iUJk,i - di [ -CI UjjUjj 



1 



(49) 



Here eij^ is the anti-symmetric Levi-Civita tensor for 
d = 3. This formula can be simplified by substitut- 
ing Tij^i — Sji — {l/2)ejimJ^m in the first term and 
'^3,1^3,1 — ^3-1^3,1 + lI'^P i'^ the second, yielding: 



fi — C£ tijkSjiLOkd — di ( -C£ SjiSji 



1 



(50) 



This is our final formula for the turbulent force. Sub- 
stituting (|50l) into (|19|l yields a similar formula for the 
circulat ion- flux : 



KeiCt) = -Ce 



CfXt) 



^ijkSjiOJkd dxi 



(51) 



According to this formula, the diffusion of vortex-lines 
out of the loop is driven by strain acting upon the gra- 
dient of the vorticity vector. This is plausible, since the 
turbulent force should act to smooth out inhomogeneities 
in the large-scale vorticity field and become negligible 
when the latter is constant. 

The same result (|51ll for the circulation-flux can be 
obtained from the "nonlinear model" of the turbulent 
vortex-force: 

/* = C£^eijkUj^iidk,i 

= Cee,jkS,iU3k,i + Qc£2|cJ|2^ . (52) 

Although this derivation yields the same result, it is the- 
oretically less well-founded because of the poorer UV- 
locality properties of the vortex-force. On the other 
hand, it gives a little more physical insight, especially 
through the following alternative expression for the 
vortex-force: 



f 



r = crv-(sxuj) + -ct{u}-v)uj 



(53) 



Here (SxaJ)ji = tikiSjk^i defines what was termed in 
[23 | the "skew-strain matrix" for d = 3. Formula H53|l 
is straightfowardly derived by calculating the divergence 
(SxU)jij and gathering the terms. This expression 
makes a nice connection with the MSG expansion for 
the turbulent stress, developed in |23|. The first term 
on the righthand side of (|53|l corresponds to one of the 
stress contributions in the MSG expansion, proportional 



to "skew-strain" . That term makes no strongly UV-local 
contribution to energy flux but a major contribution to 
helicity flux and here we see also to circulation flux. 

The second term on the righthand side of H53|) corre- 
sponds to another term from the MSG expansion in |24| . 
a contractile stress along vortex-lines, t™'^*°'^ (x —uOiUjj. 
As discussed in ^24| , the effects of the small-scale turbu- 
lence give the large-scale vortex-lines "elastic" properties. 
The second term in (|53|) therefore has a simple geometric 
interpretation and can be written as 



ic^^(aJ-V)u; 



-KCe\u3\^n. (54) 



To derive we have used the Frenet-Serret equations 
(e.g. see Q) with t = a) the unit tangent vector along 
large-scale vortex lines, n the unit normal vector and 
b the binormal |2l|. The term f™rtcx -^^ ^ 
lei to vortex lines arises from variations in the vortex- 
strcngth along the line. The term fjo^cx^ which arises 
from bending of vortex lines, is proportional to the cur- 
vature K of the line and is directed along the normal n. 
Note that ^ gives a contribution to vorticity transport, 
^vortex ^ (l/2)KC£2|uj|2b, which is directed along the 
binormal, reminiscent of the velocity of a slender vortex 
filament in the local- induction approximation . 

The formulas H50|l and H52() for the turbulent force sim- 
plify in space dimension d = 2. In that case, 



f^ 



'CfS^j{djUj)-d^{Cfa^) 
'1 



(55) 



where icr are the eigenvalues of the symmetric, trace- 
less strain matrix Sij and Sij — Sitekj = —e ikS ki is 
another symmetric, traceless matrix, called in |6l| the 
"skew-strain matrix" for d = 2. (Note that eij is the 
d — 2 anti-symmetric Levi-Civita tensor.) The corre- 
sponding result for the circulation-flux is 



Ki{c\t) = ce 



Ce{t) 



Sij{djUj) dxi. 



(56) 



This result can be derived as well from equation (|22|l and 
the "nonlinear model" for the vorticity transport vector 
in d = 2, 



ceil. 



(57) 



previously considered in j37l l38l | . (This formula is equiv- 
alent to that for the vortex-force in eg. 1)55(1 .1 Note, 
however, that the formula (|57|) predicts "down-gradient" 
transport of vorticity whenever there is a positive rate 
of vorticity-gradient stretching and this is expected in 
d = 2 both for the direct enstrophy cascade [33, HI] and 
also the inverse energy cascade |6l|,|63. "Down-gradient" 
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vorticity transport is qualitatively correct in the enstro- 
phy cascade and there H56|l may yield a good approxima- 
tion. However, in the inverse energy cascade the vortic- 
ity transport must be "up-gradient" or "anti-diffusive." 
Therefore, (|56|1 is not likely to be a good approxima- 
tion in the inverse cascade range. It must be corrected 
by higher-order terms in the convergent MSG expansion, 
corresponding to smaller subgrid scales or higher-order 
gradients. 

III. CONCLUSIONS 

The main purpose of this paper was to elaborate a 
physical theory of "circulation cascade" in classical fluid 
turbulence. We have attempted to explain the concep- 
tual basis of the phenomenon, its physical mechanisms, 
the scale-locality properties of the cascade, and its rela- 
tion to inertial-range vortex-stretching and energy trans- 
fer. Clearly, there are many important issues that call 
for further work. Chief among these is to determine the 
validity of G. I. Taylor's proposed mechanism for turbu- 
lent energy-dissipation, based on vortex line-stretching 
[19I I20I w\ . Even after seventy years of research, ba- 
sic elements of Taylor's proposal remain open to ques- 
tion. In particular, the strong inertial-range violations 
of the Kelvin Theorem — predicted in and observed 
in ^3 — cast some doubt on a key piece of Taylor's the- 
ory. It is possible that circulation conservation remains 
valid in some weaker sense, e.g. the conditional-mean 
version of the "martingale conjecture" in Further 
research is necessary to see whether any weaker form 
of the Kelvin Theorem holds at high Reynolds numbers 
and, if so, whether it is sufficient for the purposes of Tay- 
lor's mechanism. It should be emphasized that even the 
existence of circulations in the infinite-Reynolds-number 
limit is an open question. Advected loops in a turbulent 
flow are expected to become fractal 03 and defin- 
ing line- integrals for non-rectifiable curves demands some 
mathematical sophistication 17]. In superfluids the ad- 
vected contours in Kelvin's Theorem can also become 
highly distorted, with interesting consequences for vor- 
tex motion [s^. Fractality of the advected loops could 
have significant implications for conservation of circula- 
tions in fluid turbulence. 

In addition to hydrodynamics of incompressible fluids, 
there are other turbulent systems for which phenomena 



similar to "circulation cascade" are expected to exist. Of 
these, one of the most significant is magnetohydrody- 
namic (MHD) turbulence of plasmas. In this case it is 
Alfvcn's Theorem |6^ on conservation of magnetic flux 
in the ideal, zero-resistivity limit which plays the role 
of Kelvin's Theorem. However, there is strong evidence 
from observations of magnetic flux reconnection rates 
in astrophysical settings to believe that Alfven's Theo- 
rem breaks down in MHD turbulence even with negligi- 
ble resistivity |6^ . This violation of conservation of 
magnetic flux, presumably due to a similar cascade phe- 
nomenon as for Navier-Stokes dynamics, is discussed in 
a following paper for the MHD equations ■ Another 
important problem is turbulence in geophysical fluids, 
where Ertel's Theorem "68^ on conservation of potential 
vorticity (PV) plays a fundamental role in theories of 
quasi-geostrophy. It is well-known that Ertel's Theorem 
is a differential form of the Kelvin Theorem (e.g. see 
|69j . Section 2.5 or ,27,], Section 4). The "cascade of cir- 
culations" in this context should be quite similar, gener- 
ally speaking, to that for two-dimensional Navier-Stokes 
and correspond to a turbulent transport of PV out of 
the advected loop. However, in geophysical fluid dynam- 
ics there is an additional complication that the loop in 
Kelvin's theorem must lie in a surface of constant density 
(bouyancy) or pressure "G^. Thus, turbulent mixing 
of isopycnal surfaces is an additional source of breakdown 
of Kelvin's Theorem in this context. Finally, another in- 
teresting setting for "circulation cascade" is superfluid 
turbulence 70]. The analogy between quantum phase- 
slip and circulation-cascade could prove useful here. 
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